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ABSTRACT 

Eukaryotic High-Mobility Group B (HMGB) proteins 
alter DNA elasticity while facilitating transcription, 
replication and DNA repair. We developed a new 
single-molecule method to probe non-specific 
DNA interactions for two HMGB homologs: the 
human HMGB2 box A domain and yeast Nhp6Ap, 
along with chimeric mutants replacing neutral 
N-terminal residues of the HMGB2 protein with 
cationic sequences from Nhp6Ap. Surprisingly, 
HMGB proteins constrain DNA winding, and this tor- 
sional constraint is released over short timescales. 
These measurements reveal the microscopic dis- 
sociation rates of HMGB from DNA. Separate micro- 
scopic and macroscopic (or local and non-local) 
unbinding rates have been previously proposed, 
but never independently observed. Microscopic dis- 
sociation rates for the chimeric mutants (~10s~ 1 ) 
are higher than those observed for wild-type 
proteins (~0.1-1.0s -1 ), reflecting their reduced 
ability to bend DNA through short-range inter- 
actions, despite their increased DNA-binding 
affinity. Therefore, transient local HMGB-DNA 
contacts dominate the DNA-bending mechanism 
used by these important architectural proteins to 
increase DNA flexibility. 

INTRODUCTION 

High-mobility group (HMG) proteins are abundant 
nuclear proteins that support a range of cellular functions, 
from nuclear re-organization and DNA repair to cellular 
signaling (1-4). Within the hierarchy of HMG proteins 
and binding motifs, the HMGB subgroup is characterized 
by non-sequence-specific DNA binding. In eukaryotes, 



roughly 1 HMGB protein exists for every 10 nucleosomes 
(5). While histones organize DNA into chromatin, 
HMGB proteins induce bending of bare DNA, increase 
apparent DNA flexibility and have been proposed to 
weaken nucleosome structure, a necessary prelude to tran- 
scription (1,6,7). Specifically, while linker histone HI 
inhibits nucleosome sliding and promotes chromatin con- 
densation, HMGB proteins may dislocate DNA-histone 
contacts, although it is unknown whether binding is pref- 
erentially to the linker sequence or directly to the nucleo- 
some (8,9). HMGB proteins bind to the DNA minor 
groove strongly, unwinding and bending the DNA and 
enhancing the local flexibility of the double helix. 
HMGB proteins may influence transcriptional regulation 
by enhancing DNA flexibility, stabilizing looping and 
serving as chaperones for transcription factors (10-12). 

HMGB proteins carry one or two copies of a highly 
conserved 'HMG box' motif of ~80 amino acids, 
together with charged N- or C-terminal tail domains 
(13). As shown in Figure 1 [box A of human HMGB2 
and the single-box protein Nhp6Ap, an HMGB1 homo- 
logue found in Saccharomyces cerevisiae (14-17)], a single 
HMG box comprises three a helices. The L-shaped struc- 
ture is partially stabilized by hydrophobic contacts 
between the a helices, although the final structure is 
stable only on binding to DNA. These proteins bind 
into the minor groove of DNA with little sequence speci- 
ficity (10). At the binding interface lie a primary 
intercalating residue and a more weakly intercalating sec- 
ondary residue that perturb base pair stacking in the 
minor groove (Figure 1). These intercalating residues 
also endow HMGB proteins with a noted binding prefer- 
ence for bent DNA, including DNA kinked by cisplatin 
cross-linking (11,15,18) or cruciforms (19). HMGB 
cationic residues on the flexible N-terminal domain 
(NTD) or C-terminal domain often further enhance 
DNA bending, suggesting a more subtle relationship 
between local interactions and HMGB-induced flexibility 
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Figure 1. HMGB proteins studied in DNA-protein binding assays. (A) Box A of human HMGB2 bound to cisplatin-modified DNA (PDB: lctk) 
(15). Box B and the linker sequence of the tandem box structure are omitted for clarity. (B) Single-box domain of HMGB1 homologue Nhp6Ap 
from S. cerevisiae bound to a recognition sequence for the transcription factor SRY (PDB: lj5n) (14). Dark spheres illustrate partial intercalators 
methionine in Nhp6Ap, and alanine in HMGB2 (box A), as well as the intercalating residue phenylalanine present in both. (C) Sequences of 
HMGB2 (box A), Nhp6Ap and two chimeras (termed Ml and M2) based on the box A motif (residues in blue), including cationic residues from the 
N-terminus of the Nhp6Ap protein (in red — the region of substituted residues is also highlighted) (5). Shaded boxes indicate helical regions, and 
intercalating residues are shown in black. Ml and M2 were purified with poly(His) tags as shown (in grey italics), which Nhp6Ap proteins lack. Box 
A of HMGB2 was studied with and without tags. 



(20,21). HMGB-induced bending angles have been 
characterized, and estimates range from 70 to 120 degrees 
for various proteins (14,17,20-25). Enhanced apparent 
DNA flexibility may result from an ensemble of random 
hinges produced in DNA by transiently bound HMGB 
proteins (26). Atomic force microscopy analysis of the 
bending angles induced by HMGB2 (box A) showed a 
wide distribution of angles, consistent with moderately 
flexible transient hinges, rather than constant bending 
angles (27). More recent studies confirm different bend 
angle distributions for HMGB2 (box A) and Nhp6Ap 
(28). The observation of rapid association and dissociation 
kinetics in bulk experiments (26) seems to contradict 
previous observations that HMGB2 (box A) did (29) and 
Nhp6Ap did not (30) dissociate from DNA when the 
protein was removed from the solution surrounding the 
DNA molecule, suggesting slow dissociation kinetics. 

Thus, despite significant progress in our understanding 
of HMGB-DNA complexes, it is crucial to clarify the 
kinetics of these interactions. There is some evidence 
that HMGB proteins polymerize on double-stranded 
DNA (dsDNA) at high concentrations, gradually 
forming filaments over long incubation times of several 
hours (29). Other proteins have been noted to rigidify 
DNA (31,32). However, it is unlikely that this mode 
causes DNA bending and flexibility changes, which 
occur at much lower protein concentrations and on 
shorter timescales (26,29). A recent study (30) reported 



slow (10~ 2 to 10~ 3 s -1 ) DNA dissociation rates for 
several non-sequence-specific DNA-bending proteins 
(HU, Fis and Nhp6Ap) in protein-free solutions after de- 
position at low protein concentrations when no cooperative 
binding is observed. Remarkably, rates of exchange 
between bound and free proteins were found to be propor- 
tional to the free protein concentration (30). Kinetic regimes 
for HMGB/DNA interactions thus remain mysterious. 

Minor groove intercalation and cationic tail positioning 
in the opposing major groove are known mechanisms of 
DNA bending by HMGB proteins (20,21,33). However, 
the relative importance of these two contributions to 
DNA bending varies between different proteins. Here we 
compare two distinct HMGB proteins: the yeast Nhp6A 
protein with strongly cationic NTD (+7 net positive 
charge) and strongly intercalating methionine residue 
versus the human HMGB2 (box A) protein, which lacks 
a positive NTD and has a much weaker intercalating 
alanine residue (Figure 1). Intact HMGB2 is a two-box 
protein; box A binds DNA more weakly than box 
B. Nhp6Ap more closely resembles the latter box. To 
study the effects of N-terminal cationic residues, two add- 
itional constructs were prepared that graft cationic 
N-terminal residues found in Nhp6Ap in place of the 
N-terminal amino acids of the more weakly binding 
HMGB2 (box A) (5). Mutant Ml adds three positive 
charges to the N-terminus of HMGB2 (box A), whereas 
M2 adds six (Figure 1C). To characterize the roles of the 
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N-terminal charge and intercalating residue in DNA- 
binding affinity, kinetics and DNA bending, these 
proteins were exposed to single duplex DNA molecules 
held and stretched within an optical tweezers apparatus 
(34,35). We now extend our previous studies of this system 
by characterizing binding affinity for relaxed and over- 
stretched DNA, intercalation, extent of DNA bending 
and exchange kinetics. These results answer key questions 
and provide important new insights into the molecular 
mechanism by which HMGB architectural proteins 
increase the apparent flexibility of DNA while resolving 
previous seemingly contradictory kinetics measurements. 



MATERIALS AND METHODS 

Protein expression and purification 

Recombinant HMGB proteins were expressed and 
purified as described previously (36). Where the 
N-terminal hexahistidine affinity tag was removed, an 
immobilized thrombin reagent was used as described by 
the manufacturer (Thrombin CleanCleave, Sigma). 
Briefly, 200 ul of immobilized thrombin was used per milli- 
gram of fusion protein in 1-ml cleavage reactions 
incubated at 24°C for 16 h. Protein was then further 
purified by size exclusion chromatography in phosphate- 
buffered saline on a Superdex 200 10/30 column eluted at 
a flow rate of 0.4 ml/min. Desired fractions were pooled, 
dialyzed against 20 mM HEPES (pH 7.5), 100 mM KC1, 
1 mM ethylenediaminetetraacetic acid and 1 mM DTT 
and concentrated. Proteins were stored at — 20°C in this 
dialysis buffer supplemented with 50% (v/v) glycerol. 

Optical tweezers 

Phage-A. DNA, biotinylated on the opposing cohesive 
ends, is tethered between a pair of 5-um-diameter 
streptavidin-coated polystyrene beads within an optical 
tweezers flow cell, as shown in Figure 2A (34,35). 
Experimental buffers consist of 10 mM HEPES, pH 7.5, 
and lOOmM Na + . In position compliant experiments, a 
fast-scan piezoelectric stage moves the flow cell relative 
to the trap in fixed 100-nm increments. Force calibration 
is determined from the known overstretching transition 
force of DNA at specific solution and pulling rate condi- 
tions (37). Typical sets of extension and recovery data are 
shown in Figure 2C. Protein is introduced into solution 
(Figure 2B) while the tension in the DNA molecule is fixed 
at ~20 pN, to reduce the probability of protein-mediated 
loop formation (23). Cycles of extension and release are 
repeated in the presence of a given protein concentration. 
Representative cycles are shown in Figure 2D, for 400 nM 
of the protein Nhp6Ap. 



RESULTS 

Effect of HMG proteins on dsDNA elasticity 

DNA force-extension curves change significantly as 
HMGB protein concentration is increased. Figure 2 illus- 
trates the differences between free DNA and DNA in 
saturating concentrations of Nhp6Ap. Before the onset 



of the DNA overstretching transition at ~62 pN, the 
DNA and HMGB-DNA stretching curves are highly 
reproducible and near equilibrium at all pulling rates. 
Elasticity is well described by the worm-like chain 
(WLC) model. In this model, the polymer conformation 
is characterized by its contour length, B ds , persistence 
length, P ds , and optional inclusion of a phenomenological 
enthalpic stiffness (S ds ) (38,39); 



b(F) = 5, 



tlx 



1 



PdsF 



1/2 



F 

Sds 



(1) 



Fitting DNA force-extension data to the WLC model was 
first performed in the absence of protein (Figure 3A). The 
non-linear algorithm of Levenberg-Marquardt was used 
to minimize xl (40,41). Fitted values are typical and in 
reasonable agreement with previously published results 
under these solution conditions (37,38). 

Figure 3A also compares the fitted protein-free DNA 
force-extension curve with that obtained in the presence of 
saturating concentrations of Nhp6Ap. Interestingly, 
Nhp6Ap-saturated DNA is shorter than protein-free 
DNA at low forces and longer than protein-free DNA 
at high forces (>10pN). We interpret this effect as 
evidence of HMGB-induced DNA compaction by 
bending at low force (entropic elasticity regime), with 
increase of the HMGB-DNA contour length due to 
HMGB intercalation observed when the complex is 
straightened at higher forces. These effects dramatically 
change the fitted elastic parameters for DNA saturated 
by Nhp6Ap (Figure 3A). Similar results were previously 
obtained for DNA saturated a by single-box motif 
HMGB2 protein (23,29). In particular, there is a remark- 
able 6- to 8-fold reduction in the dsDNA persistence 
length, an ~10% increase in the DNA contour length 
and a minor change in elastic modulus on saturated 
binding of all HMGB proteins studied in this work. An 
alternative fitting scheme describes HMGB-DNA regions 
with the freely jointed chain model, assuming the HMGB- 
DNA complex is a polymer that contains a random 
mixture of two polymers of different elasticity: 
protein-free dsDNA and protein-saturated dsDNA, in 
Supplementary Figure SI and Supplementary Table SI. 

Interestingly, there is no evidence that these HMGB 
proteins dissociate from DNA owing to applied force, 
even up to 200 pN. Overstretching curves are reproducible 
even after several cycles of extension and release 
(Figure 2), even up to forces of 200 pN (see Figure 5). 
The ability of the applied force to straighten the 
HMG-bent DNA without protein dissociation indirectly 
supports one of the major conclusions of this work: 
HMGB-type proteins bend DNA by introducing flexible 
transient hinges, rather than rigid fixed DNA bends. This 
may reflect the ability of HMGB proteins to both compact 
DNA by bending and elongate DNA by intercalation. 
While applied force tends to dissociate proteins that 
shorten DNA (42), force promotes DNA binding for 
proteins that elongate DNA (43). Given that our data 
are reproducible over several cycles of stretch and 
release, and that changes in the measured properties 
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Figure 2. Cycles of extension and release for DNA in an optical tweezers experiment. (A) DNA is tethered within a flow cell and extended between 
functionalized spheres. (B) DNA is extended in a solution containing HMGB protein. (C) Phage-A DNA (48 500 bp in length) is extended (solid line) 
and overstretched at ~60 pN. Release (dotted line) indicates that the process is reversible on the timescales of these experiments. Three sequential 
cycles of extension and recovery (shown in red, green and blue) illustrate the reproducibility of the data. (D) When DNA is exposed to a solution 
including 400 nM Nhp6Ap, three extension and release cycles (blue, green and red) reveal consistent changes in the persistence length, contour length 
and the overstretching force. 



appear to vary smoothly with protein concentration (see 
below), it is reasonable to assume that the fractional 
HMGB saturation of DNA, 0, is approximately inde- 
pendent of force and determined by the bulk protein 
concentration. 

Figure 3B and C present the concentration dependences 
of the fitted persistence length, Pd s (c), and contour length, 
B ds {c). A simple model assumes that DNA-bending events 
are statistically independent, and independent of DNA- 
bending modes observed in the absence of protein. The 
reciprocals of the protein-free P D and protein-saturated 
P L values may be averaged, weighting by the fractional 
DNA site occupancy, 0, to yield the observed persistence 
length (P) (44,45): 



PlPd 
Pl+& ■ (Pd - Pi) 



(2) 



Variation of DNA occupancy by protein, 0, as a function 
of c is described by a binding isotherm. We applied the 
McGhee-von Hippel isotherm (46,47), conventionally 



used for the description of polymeric DNA-protein 
binding: 
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Here K D , n and co are the equilibrium protein-DNA dis- 
sociation constant, binding site size in base pairs and 
binding cooperativity parameter, respectively. The DNA- 
binding site size (« = 7) for all HMGB single-box proteins 
was estimated from crystal structures of protein-DNA 
complexes, such as HMGB (box A) presented in 
Figure 1, and others (48). The value of the HMGB DNA 
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Figure 3. Quantifying DNA-protein elasticity. (A) The extension data (solid symbols) of a typical single extension cycle for DNA (blue diamonds) 
and DNA in the presence of 400 nM Nhp6Ap (red circles). Vertical bars represent standard error of measurement when larger than symbols. Fits to 
the WLC model for bare DNA (blue line) provided values of P ds = 48 ± 2nm, B ds = 0.340 ± 0.002 nm/bp and S ds = 1400 ± 100 pN (x; = 0.45). 
When the DNA was exposed to 400 nM Nhp6Ap (red line), these fitted values became P^ = 6.7 ± 0.7 nm, B& = 0.380 ± 0.002 nm/bp and 
Sds = 1240 ± 200 pN (xl = 0.37). Residuals to these fits are shown on the right. (B) Fitted HMGB-DNA persistence length, P ds (c), determined 
versus concentration of Nhp6Ap (red), constructs Ml (green) and M2 (yellow) and HMGB2 (box A) proteins with and without an N-terminal 
poly(His) tag (blue and cyan). Symbols represent averages of fitted parameters from three to six fitted data sets. Lines indicate fits to Equations 2 and 
3, to determine the equilibrium constant, K D , and the persistence length of DNA saturated with protein, P L . (C) The HMGB-DNA contour length, 
Bd.£ c ), fit to Equations 2 and 4, to determine the contour length of DNA saturated with protein, B L . The value of K D is indistinct from the fits of 
Figure 3B. All fits assume a fixed binding site size (n = 7) and weak cooperativity (m = 20), while x 1 ,, ranges from 0.5 to 2.0. Fitted parameters for all 
proteins are summarized in Table 1. 



binding cooperativity parameter that best fits P,is( c ) an d 
B ds (c) dependencies was found to be co = 20 for all fits, 
corresponding to a free energy of protein-protein inter- 
action, ksT-lnco, of 1.7kcal/mol. This indicates weak 
cooperativity of HMGB binding, as previously observed 
(27,49). The persistence length of bare DNA (P D ) was 
fixed at 45 nm, such that the only free fitting parameters 
were P L and K D . Best-fit values of these parameters are 
then summarized in Table 1. An alternative fitting 



procedure using the Hill binding isotherm, presented in 
Supplementary Figure S2 and Table S2, yields Hill par- 
ameters ranging from 1 to 1.5, also consistent with weakly 
cooperative HMGB-DNA binding. Variations in any 
results due to the presence of the his 6 tags were minor, 
and these differences are considered explicitly in 
Supplementary Tables S3 and S4. 

The DNA contour length increases by up to 10% in the 
presence of HMGB proteins, confirming the intercalative 
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nature of HMGB-DNA binding. While this is a relatively 
small effect compared with the major reduction of persist- 
ence length by HMGB proteins, it implies DNA elong- 
ation by the equivalent of a single base pair (0.34 nm) 
for every 10 bp. Taking into account that each HMGB 
protein binds to at least 7 bp and has a single primary 
intercalating aromatic residue, our result suggests an 
almost complete intercalation of this residue, comparable 
with the DNA elongation on intercalation of conventional 
intercalators, such as ethidium and intercalating ruthe- 
nium complexes (43). Quantitatively, the DNA contour 
length increases owing to independent HMGB-binding 
events that add to the overall length, in proportion to 
their occupancy (0): 

Bdsiy) = Bd+& ■ (B L - B D ) (4) 

As above, the binding site size is fixed at n = 7 bp, 
co = 20 and B D = 0.34nm/bp. This leaves only K D and 
B L as free parameters for the fits to Equations 2 and 4, 
which are shown in Figure 3C, and collected in Table 1. 

Using models of polymer elasticity, the decrease in the 
apparent DNA persistence length due to the protein 
binding can be further related to the average DNA- 
bending angle, p, associated with individual protein- 
binding events. Assuming independent binding/bending 
events and taking into account that one protein binds 
per ri'B L length of dsDNA at saturation, p can be 
related to the measured saturated persistence length 
P L (44,45,50,51): 

Calculated p values are summarized in Table 1. 

The fact that Kj values obtained by fitting two inde- 
pendently measured characteristics of the HMGB-DNA 
complex, Pd s (c) and B ds (c), are in excellent agreement 
confirms our ability to characterize protein binding 
under applied force. Fitted binding affinities for all four 
HMGB proteins summarized in Table 1 correlate with the 
positive charge of the NTD as previously suggested (5). In 
contrast, saturated values of the HMGB-DNA complex 
persistence length P L , contour length B L and bending 
angle per protein p do not correlate with NTD positive 
charge (Table 1). Similar 10-fold reductions in DNA per- 
sistence length are induced by the two wild-type proteins, 
HMGB2 (box A) and Nhp6A, despite different NTD 
charges. The effects of the two chimeric HMGB proteins 
with altered NTD residues are approximately a factor of 
two weaker. 

We previously observed the formation of a rigid 
HMGB-DNA filament at much higher HMGB concen- 
trations (>10uM) and on longer timescales of several 
hours (29). Filament formation has also been observed 
for HU, a functional analog of HMGB found in 
bacteria (31,52). However, HU-induced dsDNA bending 
and compaction were directly followed by dsDNA 
straightening as HU concentration was increased. 
Filament formation at HU concentrations leading to 
DNA straightening resulted in slower HU dissociation 
kinetics (52). In contrast, as shown in this work, increasing 



HMGB concentration well below 10 uM leads to saturated 
protein binding and a 10-fold reduction in DNA persist- 
ence length without filament formation. Thus, HMGB 
effects on DNA discussed here are not directly related to 
HMGB-HMGB interactions inducing filament formation 
at much higher HMGB concentrations. 

Effect of HMGB proteins on equilibrium DNA 
strand pairing stability 

Thus far, we have considered HMGB effects on DNA 
elasticity before the onset of the force-induced overstretch- 
ing transition in dsDNA. During this transition, the two 
DNA strands are unwound and stretched to almost twice 
their B-form length. The nature of the overstretching tran- 
sition is still under investigation (53-58), but the observed 
increase in overstretching force required in the presence of 
protein binding indicates preferential HMGB stabilization 
of B-form dsDNA (Figure 4A). The DNA overstretching 
transition appears to be reversible in the presence of 
HMGB protein, such that DNA stretching and release 
curves overlap, and are independent of the pulling rate 
when the DNA is stretched at the rate of lOOnm/s or 
slower. Because the total X-DNA extension increases by 
~10 urn during overstretching, equilibrium is preserved on 
the timescale of ~100 s or longer. We now analyze HMGB 
protein effects in this overstretching regime to gain 
important new insights. 

For this purpose, we define the overstretching force, F ov , 
as a transition midpoint force averaged between DNA 
extension per base pair values of 0.42 and 0.48 nm. In 
the absence of proteins, F or is known to vary with 
solution ionic strength, pH and temperature (37,59,60). 
Under our experimental conditions (100 mM Na + , 21°C, 
pH 7.5), F m , = 62.6 ± 0.4 pN, in accord with previous 
studies. Figure 4A shows that a saturating Nhp6Ap con- 
centration increases F av to ~73pN. This increase in 
the transition force characterizes the equilibrium preferen- 
tial binding of HMGB protein (per base pair) to 
B-form DNA, relative to the overstretched state: 
SG = 8F- Ax % 0.5 kgT. For a binding site size of n = 7, 
this implies preferential HMGB protein binding to B-form 
DNA that is a factor of e " A * « r / k i> T % e 3 » 20 stronger than 
to overstretched DNA. This factor of 20 preference for 
B-form DNA (which is similar to the measured 
cooperativity parameter co) may simply reflect the loss of 
cooperative HMGB-HMGB binding interactions on 
overstretching without the release of HMGB-DNA 
contacts. Complete DNA overstretching is typically 
observed in the presence of bound protein which argues 
against filament formation (Figure 2D). The high repro- 
ducibility of DNA stretch-release curves even under 
saturating protein concentrations suggests that HMGB 
proteins remain bound throughout the stretching cycle, 
accommodating the transition from B-form DNA to over- 
stretched DNA on a timescale longer than 100 s. 
Additional support for this hypothesis comes from the 
fact that faster DNA overstretching reveals non- 
equilibrium protein effects on the overstretching transi- 
tion, as discussed below. It is exactly this insight that 
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Figure 4. Quantifying DNA-protein stability. (A) The overstretching 
force for a DNA molecule (blue circles, triangles and diamonds) rises 
when that molecule is exposed to 400 nM Nhp6Ap (red circles, triangles 
and diamonds). The extension data of three contiguous cycles of ex- 
tension (solid symbols) and release (open symbols) are plotted to illus- 
trate the reproducibility of the data. Averaging extensions over the 
range of the graph (0.42-0.48 nm per base pair), and over each of the 
three cycles, the overstretching force is found to be 62.4 ± 0.4 pN in 
the absence and 70.2 ± 0.4 pN in the presence of 400 nM Nhp6Ap. 
(B) The observed overstretching force F ov (c) fit to the site exclusion 
binding isotherm of Equations 3 and 6 to determine K D and the over- 
stretching force of DNA saturated with protein, F^ r All fits assume a 
fixed binding site size (n = 7) and weak cooperativity (co = 20), while x 2 v 
ranges from 0.5 to 2.0. Fitted parameters are found in Table 2. 



permits estimation of DNA dissociation times for 
individual HMGB proteins. 

Assuming that F ov (c) increases linearly with the amount 
of HMGB protein bound to DNA, the complete F ur (c) 
dependence presented in Figure 4B can be fit to the 



binding isotherm given by Equation (2) using the 
expression 

F m (&) = FZ+& ■ (F L ov - F D m ) (6) 

Here F^ r and F% are the protein-free and protein-saturated 
values of F ov . We fit the experimental F ov (c) to Equation 6, 
with ©(c) given by the McGhee-von Hippel isotherm 
(Equation 3). Fits appear in Figure 4B, and the resulting 
equilibrium dissociation constants (K D ) and saturated 
melting forces (F^ v ) are shown in Table 2. Interestingly, 
the K D values obtained from F ov (c) fits are ~4-fold higher 
than the K D values obtained by fitting the P</ S (c) and B c / S (c) 
dependencies (compare K D values in Table 1 and 2). There 
are many possible reasons for this minor discrepancy, such 
as an inaccurate description of the protein effect on the 
transition force by the empirical relationship given in 
Equation 6. Fits of F m ,(c) to the Hill isotherm are shown 
in Supplementary Figure S3 and Table S5. Most import- 
antly, consistent K D variations between the four HMGB 
proteins studied in this work suggest that in all cases, we 
are measuring the characteristics of the same binding 
process. K D decreases (reflecting stronger binding affinity) 
with increasing charge on the NTD of the HMGB pro- 
tein, in accord with the previous bulk solution measure- 
ments (5). 

HMGB proteins torsionally constrain B-form DNA, 
inhibiting unwinding and re-winding on short timescales 

As discussed above, equilibrium values of F ov (c) in the 
presence of HMGB proteins can only be observed if the 
DNA stretching time is >100s, i.e. when the DNA is 
pulled at the rates slower than lOOnm/s. However, it is 
possible to probe DNA-ligand interactions on much 
shorter timescales by measuring forces at much faster 
DNA extension rates. This is accomplished by decreasing 
the signal averaging of the detector, allowing pulling rates 
to increase up to 50000nm/s, thereby reducing the total 
DNA stretching time to ~0.2 s. In the presence of HMGB 
proteins, cycles of extension and recovery at increasing 
pulling rates show increased overstretching forces 
relative to bare DNA. Figure 5A and B show responses 
for HMGB2 (box A) and the mutant Ml at 
near-saturating concentrations of protein (c > K D ). 
Additional data are shown in Supplementary Figure S4. 
As the rate of extension increases, the DNA transition 
force grows steeply from the protein-saturated equilibrium 
value Fg V ~ 73 pN to F T oy ~ 115 pN, far more than the 
transition force for bare DNA (Figure 5C). In the absence 
of protein, there is a weak dependence of the overstretch- 
ing force on the pulling rate, and DNA remains close to 
equilibrium during overstretching (61). In the presence of 
HMGB, the transition force changes over a relatively 
narrow range of pulling rates around some characteristic 
rate, v HMG , which is unique for each protein. The 
maximum transition force, F^, ~ 115 pN, is the same for 
the four HMGB proteins studied in this work and is 
similar to the overstretching force observed for bare but 
torsionally constrained DNA (55,62). 

While the stretching transition force (~ 115 pN) is well 
defined for sufficiently high protein concentration 
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Figure 5. Probing the dynamics of HMGB-DNA binding. (A) Cycles of extension and release (solid and dotted lines) for (A) 200 nM HMGB2 (box 
A) and (B) 200 nM Ml as a function of the rates shown (additional data can be found in Supplementary Figure S4). At pulling rates higher than the 
natural dissociation rate, DNA becomes torsionally constrained. Release to lower DNA extensions shows multiple plateaus, described in the text. 
(C) The average overstretching force as a function of the experimental pulling rate, v, for bare DNA (black) and DNA exposed to Nhp6Ap (red), 
constructs Ml (green) and M2 (yellow) and HMGB2 (box A) (blue) proteins. Data are averaged over active concentrations (c > K D ) as described in 
the text (and Supplementary Figure S5). Dotted lines are fits to Equation 8, giving the dissociation rate, k„jf. Rates from these fits are shown in 
Table 2. (D) Fits to Equation 1 measure P ds for bare DNA (black), 500 nM Nhp6Ap (red), 2500 nM M2 (yellow), 4000 nM Ml (green) and 200 nM 
HMGB2. Dotted lines mark values at the lowest pulling rate (100 nm/s). (E) The measured stiffness, 5, of the HMGB-DNA complexes after 
overstretching, at low pulling rates (<1000nm/s). (F) The fraction of DNA that is melted by force,/, decreases with increasing protein concentration, 
as dsDNA is progressively converted into a stable filament that cannot be overstretched at the observed forces. 



Nucleic Acids Research, 2013, Vol. 41, No. 1 175 



Table 1. Comparisons of the flexibility of DNA in the presence of HMG single-box binding proteins, arranged by net charge of the N-terminus 
(the substituted residues of Figure 1C) 

Protein q N ?i(nni) B^nm/bp) S L (pN) K D (nM) P(°) 



HMGB2 (box A) a -1 4.6 ± 0.5 0.378 ± 0.001 1600 ± 200 180 ± 30 57 ± 4 

Ml 3 11.3 ± 1.1 0.360 ± 0.002 1600 ± 200 170 ± 30 37 ± 4 

M2 7 10.0 ± 1.0 0.359 ± 0.002 1200 ± 200 81 ± 16 39 ± 4 

Nhp6Ap 7 5.5 ± 0.5 0.378 ± 0.002 1000 ± 200 71 ± 14 52 ± 5 



The values of the persistence length, contour length and stiffness of DNA in saturating concentrations of ligand are shown as P L , B L and Sl, 
respectively. These parameters and the equilibrium dissociation constant (K D ) were determined by fitting the concentration dependence of the 
persistence length, Pj s ( c ), an d contour length, B ds (c), to Equations 2, 3 and 4, with fits shown in Figure 3B and C (the results for K D are the 
same within uncertainty). The binding site size was assumed constant for all proteins at n = 7, and the cooperativity was also fixed at to = 20. Mean 
bending angles (P) were calculated from Equation 5. 

*K D and p are values averaged for HMGB samples prepared with and without poly(His) tags, which showed only minor differences. Values for 
proteins with and without the tags are shown distinctly in Supplementary Table S3. 



Table 2. Comparisons of the stability and kinetics of HMGB-DNA interactions, arranged by net charge of the N-terminus (q N ) 



Protein q N ^ v (pN) jf„(nM) k o/) (s^') /t a (xl0 8 M _1 s _1 ) 



HMGB2 (box A) a -1 73.3 ± 0.2 880 ± 70 0.88 ± 0.09 0.010 ± 0.002 

Ml 3 72.3 ± 0.3 480 ± 60 3.3 ± 0.5 0.069 ± 0.017 

M2 7 72.2 ± 0.3 260 ± 30 6.1 ± 0.7 0.23 ± 0.05 

Nhp6Ap 7 73.3 ± 0.3 160 ± 20 0.15 ± 0.01 0.009 ± 0.002 



Equilibrium dissociation constants (K D ) and HMGB saturated overstretching forces (F^ v ) are determined from the change in F ov (c), fitted to Equation 
3 and Equation 6 as illustrated in Figure 4B. Rate dependence experiments resolved HMGB-DNA unbinding rates (k of f) using fits to Equation 8, 
shown in Figure 5C. Unbinding rates and association rates (k a = k 0 jf/Ko) are deduced for each protein. 

ll K D values averaged for HMGB prepared and without poly(His) tags. Values for proteins with and without the tags are shown distinctly 
in Supplementary Table S4. 



(c > K D ), DNA release curves vary. Sometimes the DNA 
release force plateau coincides with the stretching plateau 
at 115pN. In other cases, however, the strand 
re-annealing transition occurs at a much lower (but 
well-defined) force of ~50 pN that is independent of the 
HMGB protein type (Figure 5A and B, see also 
Supplementary Figure S4). Furthermore, some DNA 
release curves begin with the strand re-annealing plateau 
at the higher 1 1 5 pN force, but then jump to the lower 
force of 50 pN. These double-plateau DNA stretching and 
release curves were previously observed during the 
stretching of torsionally constrained bare dsDNA where 
the degree of unwinding was fixed by immobilizing both 
DNA strands at both termini. DNA unwinding was then 
achieved by rotating the end-fixed dsDNA around its axis 
by a known number of turns before stretching in a 
magnetic tweezers apparatus (62). Relative lengths of the 
low and high force plateaus were determined from frac- 
tional DNA unwinding, with the transition for torsionally 
constrained B-form DNA at 1 15 pN and the transition for 
totally unwound dsDNA overstretching at 50 pN. Non- 
equilibrium DNA under-winding often persists until 
complete DNA release, leading to significant hysteresis in 
the DNA stretch-release cycle up to small forces below 
50 pN. Based on HMGB-saturated dsDNA force- 
extension curves, we hypothesize that HMGB proteins 
inhibit DNA unwinding and re-winding at high pulling 
rates. When v >> v HMG , both the extension and relaxation 
curves plateau at 115 pN with the DNA fully constrained. 
At rates near the characteristic rate of the protein, 
v ~ v hmg, some slipping and release of strand winding 
occurs after complete DNA overstretching and on the 



onset of DNA release. This leads to the observed double 
plateau in the DNA release portion of the cycle, during 
which the DNA is trapped in a partially underwound 
state. In other cases, when such winding relaxation at 
maximum DNA extension is complete, the single DNA 
release plateau at ~50 pN is observed. In cases of incom- 
plete DNA saturation with HMGB protein (at c < K D ), 
decreasing the pulling rate below v HM g leads to a 
decrease in the force of the extension plateau, whereas the 
force of the release plateau rises, with both forces finally 
converging at the equilibrium value of ~73 pN at low 
pulling rates, v << v HMGB . This observation suggests a 
partial HMGB protein-induced constraint on DNA un- 
winding and re-winding by protein concentrations below 
saturation. On stretching such partially torsionally con- 
strained DNA, some unwinding may occur. However, 
this rate is lower than the rate of DNA unwinding 
imposed by overstretching, leading to increased transition 
force with pulling rate. These considerations suggest that 
on short DNA pulling timescales (0.2-10 s), HMGB 
proteins do not interfere with the DNA overstretching 
transition, but create a block to unwinding and re- 
winding of DNA as it is overstretched. 

Importantly, the pulling rate dependence of the DNA 
overstretching force is only observed at higher protein 
concentrations (c>K D ) where the behavior of F ov (v) 
becomes highly reproducible and independent of the 
history of the particular DNA molecule. This observation 
rules out enhanced protein binding or un-binding at high 
forces. The effects of protein concentration are illustrated 
further in Supplementary Figure S5. DNA winding is 
apparently restricted by individual proteins bound to 
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both B- and overstretched forms of DNA. Complete tor- 
sional constraint would require at least one protein per 
helical turn of dsDNA (1 protein per 10 bp). Taking 
into account that each single-box HMGB-type protein 
covers 7 bp, this requirement amounts to fractional 
DNA binding of 0 >0.7, leading to the requirement of 
near-complete DNA saturation. It is important to note 
that complete DNA coverage by HMGB is not required 
for observation of the HMGB-induced DNA torsional 
constraint. In fact, the protein filament formation on 
DNA would necessarily interfere with DNA overstretch- 
ing, leading to the altered length of overstretched DNA 
and/or irreproducible overstretching behavior resulting 
from random filament breaking pattern, as discussed 
below. If the binding of individual HMGB proteins 
indeed restricts DNA strand winding, the time for 
winding relaxation will become limited by the unbinding 
time of such proteins. In the context of an optical tweezers 
experiment, this is the time required to overstretch DNA 
at a given pulling rate, v HMG = NAx/t. The change in 
length per base pair on overstretching in the presence of 
HMGB is Ax = Xds~x ss = 0.18nm/bp, while the total 
number of 1-DNA base pairs is TV = 48 500. An observ- 
able effect on DNA winding is anticipated when the 
protein dissociation time becomes similar to the DNA un- 
winding (overstretching) time. The desired protein dissoci- 
ation rate of interest ik 0 ff) can be thus found from the 
characteristic unwinding rate: 



Kff = vhmgb/NAx 



(7) 



Fitting the experimental F ov (v) dependencies for each 
protein presented in Figure 5C to a simple sigmoidal 
function determines the midpoint rate, v HM gb- 



F„ v (y) 



v+ vhmgb 



(F T - F L ) 

V ov ov/ 



(8) 



Here F% = 73 pN, and F^ r =115 pN. Combining the dis- 
sociation rate k a ff (calculated according to Equation 7 
using the fitted values of v HMGB ) with the equilibrium dis- 
sociation constants (K D ) determined from the F ov (c) fit to 
Equation (6) allows us to estimate the association rate 
(k a = k 0 ff/K[)). These values are summarized in Table 2. 
Calculated k a ff values for all four HMGB proteins studied 
are relatively fast, ranging between 0.1 and 10s _1 , where 
the two wild-type HMGB proteins are the slowest 
[0.88 ±0.09 and 0.15 ± 0.01s -1 for HMGB2 (Box A) 
and Nhp6Ap, respectively], and the two chimeric 
HMGB proteins Ml and M2 are faster (3.3 ± 0.5 and 
6.1 ± 0.7 s" 1 , respectively). Although these rates are fast, 
they are much slower than the unperturbed rate of 
protein-free DNA unwinding from free ends [~10 6 s _1 
(63)] or DNA unwinding rates limited by the rotational 
friction of A-DNA (61). We conclude that it is the actual 
dissociation kinetics of individual HMGB-DNA molecules 
that limit the DNA unwinding rate, and use this new insight 
to resolve an apparent contradiction with previously 
measured HMGB and Nhp6Ap kinetics (29,30). 

Our observation of rapid HMGB dissociation rates 
(complex lifetimes of 0.1-1 0 s) is in apparent contradiction 
with the much longer times (~ 1000 s) required to elute 



HMGB proteins from DNA by washing (42). We hy- 
pothesize that the fast dissociation rates measured in our 
DNA overstretching assays correspond to the microscopic 
dissociation rate, k 0 ff micro ~ 0.1-10 s -1 (42), in contrast to 
the macroscopic HMGB protein dissociation rate (42,64). 
Macroscopic dissociation describes the departure of 
protein from DNA into bulk solution. In contrast, micro- 
scopic dissociation involves the breaking of only 
short-range protein-DNA contacts, with the protein 
remaining within the ~l-nm range of strong electrostatic 
attraction (65). The latter process facilitates re-association 
of the same protein with DNA, while also allowing for 
protein exchange with other molecules in the bulk 
solution. We reason that microscopic HMGB-DNA dis- 
sociation events are sufficient to release the HMGB- 
induced torsional constraint on dsDNA unwinding and 
re-winding. 

The double-box HMGB2 (box A + B) also exhibits rapid 
DNA association-dissociation kinetics 

Interestingly, the double-box HMGB2 (box A + B) studied 
extensively in our previous work caused a comparable 
increase in the DNA overstretching force with increasing 
pulling rates, although it was not analyzed at the time (23). 
We have now analyzed these previous data as described 
above, and summarize the results in Supplementary 
Figure S6 and Table S6. Importantly, the double-box 
protein HMGB2 (box A + B) displays a dissociation rate 
similar to the wild-type HMGB2 (box A) and Nhp6Ap 
proteins. The association rate of HMGB2 (box A + B) is 
~10 7 M _1 s _1 , much faster than the association rate of 
HMGB2 (box A), leading to an ~ 100-fold stronger 
binding affinity for the HMGB2 (box A + B) compared 
with HMGB2 (box A) and Nhp6Ap proteins. We 
conclude that the double-box HMGB proteins exhibit 
rapid dissociation kinetics, similar to that of the single- 
box HMGB proteins, with two boxes contributing to net 
~ 100-fold stronger binding and the analogous faster 
association rate. 

DNA appears stiffer at higher pulling rates 

As the melting force in the presence of HMGB grows with 
the pulling rate, so does the apparent DNA persistence 
length. Analysis of Figure 5A and B rapid stretching 
data (before the overstretching transition) was performed 
by collecting several (typically 3) extension curves and 
averaging. The resulting values of Pd s {v) are shown for 
bare DNA and for saturating concentrations of HMGB2 
(box A), Ml, M2 and Nhp6Ap (Figure 5D).The persist- 
ence length shows no variation with the pulling rate for 
bare DNA. In the presence of HMGB protein at fast 
pulling rates (v = 50000nM/s), the fitted persistence 
lengths increase. The maximum increase in P L compared 
with its equilibrium minimal value reaches no more than 
30% even at the highest pulling rates. Apparent stiffening 
of HMGB-bound DNA on short timescales (<1 s) could 
result from unwinding inhibition imposed by bound 
HMGB on DNA or possibly insufficient time for DNA 
bending by the HMGB protein. 
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HMGB remains bound to overstretched DNA strands 
without affecting their stretching behavior 

After overstretching of HMGB-bound DNA, force 
increases steeply with extension, just as in the absence of 
the protein. The slope of this curve (S; see labels in Figure 
5A) provides a measure of the stiffness of the over- 
stretched HMGB-DNA complex. The stiffness is 
graphed for several HMGB protein concentrations in 
Figure 5E and compared with the value of 5 for bare 
DNA. For these protein concentrations and pulling 
rates, there is no measurable protein effect on the elasticity 
of overstretched DNA. This result is unexpected because 
HMGB proteins remain bound to two overstretched DNA 
strands as indicated by interference with strand unwinding 
and re-winding on timescales shorter than the dissociation 
time of the protein. Thus, HMGB binding to over- 
stretched DNA strands must be compliant, as binding 
persists after overstretching without detectably affecting 
elasticity. 

High HMGB concentrations promote formation of 
an HMGB filament that prevents overstretching 

Some DNA stretching curves display shorter overstretch- 
ing plateaus, particularly at protein concentrations 
> 1 uM. Figure 5F shows the overstretchable fraction of 
the total DNA length as a function of protein concentra- 
tion, calculated as the apparent DNA overstretching 
length per bp, Ax, divided by its maximum value, 
Bss-Bjs, <f= Ax/(B,,,-B ds ), where B ss -B ds = 0.22 nm/ 
bp). The overall magnitude of this effect is modest in the 
studied range of protein concentrations (up to 4uM), 
reaching at most 20%. 

DISCUSSION 

Rapid HMGB-DNA microscopic dissociation versus 
slow macroscopic dissociation 

A principal clarifying contribution of this work is the 
measurement of rapid DNA dissociation rates [~0.1 to 
~6s _1 (Table 2)] for all four HMGB proteins studied 
here. These rapid DNA dissociation rates agree with 
recent bulk measurements (26) but contrast with the 
slow dissociation rates (10~ 2 to 10~ 4 s -1 ) observed for 
the same protein preparations on washing (52). Slow dis- 
sociation rates were also recently observed for Nhp6Ap 
and DNA-bending proteins HU and Fis (32,52). 
Explanations for slow dissociation kinetics included 
protein-protein-DNA binding cooperativity or sequence- 
specific protein-DNA interactions. These explanations are 
not pertinent for the HMGB proteins studied here. We 
appeal rather to the existence of two distinct protein- 
DNA dissociation rates: macroscopic and microscopic. 
Macroscopic dissociation allows proteins to escape into 
bulk solution. Microscopic dissociation breaks immediate 
short-range protein-DNA contacts, but the protein 
remains 'territorially' bound to DNA through long-range 
electrostatic interaction. Support for a microscopic 
dissociation process comes from the observation that 
slow macroscopic dissociation is accelerated by high salt 



for HU, Fis and HMGB proteins (32). The existence of 
this process also explains the observation that such 
proteins can be removed from DNA within seconds by 
the addition of excess free DNA (52), or exchanged with 
other proteins in solution at a rate proportional to the 
protein concentration (30). Previous bulk experiments 
suggested that the macroscopic dissociation rate could 
be dramatically reduced simply by reducing the concen- 
tration of free HMGB protein through addition of excess 
carrier DNA (5). As noted by Marko and coworkers (32), 
there is a notable contradiction between the measured 
slow protein-DNA dissociation rate, k 0 jj macr0 M0~ 2 to 
10" 4 s" 1 , the apparent K D of 10-100 nM and the 
observed fast DNA binding at concentrations c ~ K D . 
The association time r on (0.1-10 s) suggests a bimolecular 
association rate of k a = \/{x 0 „ ■ K D ) ~10 6 -10 8 M _1 s _1 . 
The expected dissociation rate for such binding 
k 0 ff— k a ■ Kq = l/r o „~0.1-10s _1 is much faster than the 
observed dissociation rate into solution, discussed above. 
We now measure HMGB-DNA dissociation rates that are 
in good agreement with the expectations for a bimolecular 
reaction. We suggest that these rates reflect the micro- 
scopic dissociation rates, due to the loss of HMGB- 
DNA short-range interactions. Breaking such contacts 
facilitates changes in DNA strand winding. Thus, our 
novel single-molecule experimental approach measures 
the characteristic rate at which force-induced DNA 
unwinding in the presence of HMGB becomes independ- 
ent of the unwinding rate, yielding estimates of the micro- 
scopic dissociation rate for individual HMGB proteins. 

Short range HMGB-DNA contacts play the dominant 
role in HMGB-induced DNA bending 

Our measured A: offimicro values of 0.1-10 s _1 for the four 
HMGB proteins studied here span almost two orders of 
magnitude, with the two slower rates of 0.88 and 0.15 s _1 
for HMGB2 (box A) and Nhp6Ap, respectively, and the 
two faster rates of 3.3 and 6.1 s -1 for the chimeric HMGB 
proteins Ml and M2, respectively. As A: 0 ff ]in i cro does not 
reflect the total strength of protein-DNA binding, but 
rather pertains to the barrier for breaking only local inter- 
actions, it is reasonable that the value of k o{ ^ micio does not 
correlate with apparent K D , which decreases with the net 
positive charge of the NTD (Tables 1 and 2). This result 
reflects the fact that short-range HMGB-DNA contacts 
are stronger for the two wild-type proteins than the two 
chimeric proteins. Local contact strength may depend on 
protein charge in a complicated manner that depends on 
positioning particular HMGB amino acids within the 
DNA structures. 

The fc 0 ff,micro rates correlate directly with the ability of 
HMGB proteins to compact DNA, characterized in this 
study by the smaller HMGB-saturated DNA persistence 
length, and with the ability of HMGB to extend the DNA 
contour length by intercalation (Tables 1 and 2). Thus, 
DNA bending by HMGB proteins is more correlated 
with the strength of local DNA contacts than with 
N-terminal positive charge. This finding relates to the 
work of Privalov and coauthors (66-68), decomposing 
the energetics of HMGB-induced DNA bending into 
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contributions from short-range interactions versus 
long-range electrostatic effects. Here we observe that a 
sequence non-specific HMGB protein with the smallest 
positive NTD charge, HMGB (box A), induces the 
largest DNA bend, consistent with stronger local 
contacts. This result suggests that short-range non- 
electrostatic interactions are primarily responsible for the 
DNA bending observed for the series of four proteins 
studied here. 

The slower dissociation rate of the Nhp6Ap protein 
compared with HMGB2 (box A) measured in this work 
is consistent with the superior facilitation of DNA cycliza- 
tion by Nhp6A in vitro and in vivo (5). Interestingly, 
Maher and coworkers observed that the chimeric 
HMGB proteins Ml and M2 studied in this work were 
almost as effective in promoting DNA cyclization as 
Nhp6Ap, and much more effective than HMGB2 (box 
A). This result does not contradict our current conclusion 
that the M 1 and M2 proteins promote less DNA bending, 
as the DNA-binding affinities of these chimeric proteins 
are higher, which under conditions of low protein satur- 
ation may dominate the overall protein effect on 
cyclization. 

HMGB proteins bind two DNA strands of any secondary 
structure and winding with fast microscopic kinetics 

Measurement of rapid HMGB-DNA microscopic associ- 
ation and dissociation rates depends on recognition of the 
existence of two kinetic regimes in force-extension experi- 
ments. When the complete DNA overstretching and 
release cycle is performed on a timescale longer than 
100 s, the resulting HMGB-DNA force-extension curves 
appear different than for bare DNA, while displaying re- 
producibility between cycles and lacking significant hys- 
teresis (Figure 2C and D). As increasing force changes the 
DNA conformation and the HMGB-DNA binding mode, 
the effective elastic behavior of the complex, nonetheless, 
remains well equilibrated, suggesting that the HMGB- 
DNA association-dissociation kinetics are faster than 
the DNA stretching process. This situation changes as 
the time of HMGB-DNA complex stretching becomes 
<10s. Before the overstretching transition, changes in 
the HMGB-DNA force-extension curve with pulling 
rate are moderate and appear as an increase in persistence 
length at the highest pulling rates (Figure 5D). There are 
many possible reasons for this effect, including change in 
the HMGB-DNA stretching behavior due to the 
protein-imposed torsional constraint. This constraint 
becomes more evident as a saturated HMGB-DNA 
complex is pulled through the overstretching transition 
at a high rate. The presence of the DNA-bound HMGB 
protein is only apparent through the effective restriction 
on DNA unwinding on stretching and DNA re-winding 
on release. The overstretching plateaus for both torsion- 
ally constrained B-form DNA and completely unwound, 
overstretched and torsionally constrained protein-bound 
DNA agree with their protein-free values at ~ 1 1 5 and 
~50pN (62). Furthermore, the length and elasticity of 
the overstretched DNA remain largely unaffected by 
protein (Figure 5E and F). 



These observations imply that all four studied HMGB 
proteins remain bound to DNA throughout overstretch- 
ing, strongly interfering with any changes in the linking 
number of the two DNA strands on short timescale of 
0.1-10 s. This result reflects HMGB protein promiscuity 
in binding to DNA strands with different secondary struc- 
tures. Although this result was unexpected, it is consistent 
with the known ability of HMGB proteins to bind many 
deformed DNA structures, often with greater affinity than 
to B-form DNA. Indeed, HMGB proteins are known to 
bind tightly to bent DNA (15), cross-linked DNA (33,69), 
extruded cruciform DNA (19), locally mismatched and 
unpaired DNA (70) and supercoiled DNA (71). HMGB 
proteins thus endow DNA with a flexible hinge-bending 
mode (27,28) and local softening of DNA torsional 
rigidity, as was previously observed for the HU protein 
(32). According to the present work, HMGB proteins 
transiently introduce such highly deformable sites with 
microscopic dissociation rates of 0.1-10 s -1 . 

Functional implications of the ability of HMGB proteins 
to strongly and transiently deform DNA 

The observed rapid microscopic rates of HMGB-induced 
DNA conformational interconversion facilitate cycliza- 
tion of short (100-300 bp) DNA fragments and synapsis 
of DNA sites for transcription factor binding (72). Indeed, 
HMGB proteins strongly enhance these processes even at 
low protein concentrations (c « K D ), reflecting in vivo 
stoichiometrics of one HMGB protein present per several 
hundreds of genomic DNA base pairs (5). Under such low 
HMGB-binding conditions, a DNA fragment of the 
length comparable with its unperturbed persistence 
length of 50 nm (150 bp) may be bound, on average, by 
just a few proteins. Such short DNA fragments with few 
HMGB-induced bends display greatly reduced persistence 
lengths (72,73). This behavior suggests that the few 
HMGB-induced DNA bends are not static, but 
fluctuating on the short timescale of 0.1-10 s, leading to 
an effective time-averaging rather than length-averaging 
of the observed DNA-protein complex persistence length. 

The ability of HMGB proteins to strongly but transi- 
ently bind a wide variety of dsDNA structures may facili- 
tate destabilization of nucleosomes as required for 
transcription (6,72). Indeed, DNA wound onto nucleo- 
somes is strongly deformed (74,75). The promiscuous 
HMGB-binding properties discussed above may introduce 
sites of strong DNA deformability. Based on the results of 
this study, such deformable DNA sites will not be static, 
but transient on the timescale of 0.1-10s. This behavior 
will soften both structured and unstructured DNA. 

The newly discovered ability of the HMGB proteins to 
constrain DNA winding on short timescales suggests 
novel physiological functions such as helicase retardation. 
The present study suggests that an HMGB-bound turn of 
B-form DNA can unwind on a timescale not faster than 
several seconds. In the limiting case of HMGB-DNA sat- 
uration, the rate of helicase-driven duplex unwinding 
might then be decreased to 1-10 bp/s. HMGB-type 
protein binding might slow helicase-driven duplex 
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unwinding 10- to 100-fold relative to the unwinding rate of 
free DNA. 



SUMMARY 

In this work, we stretch single DNA molecules in the 
presence of four different HMGB proteins to characterize 
HMGB-DNA interactions. In contrast to previous 
micromechanical studies of other DNA-bending 
proteins, we stretch DNA completely through its over- 
stretching transition and follow HMGB effects on both 
B-form and overstretched DNA. Importantly, we 
perform a systematic study of the effect of pulling rate 
on HMGB-DNA stretching curves, allowing characteriza- 
tion of both equilibrium and kinetic parameters for 
HMGB-DNA interactions. Slow DNA stretch-release 
cycles (~100s) appear to be in equilibrium, suggesting 
that HMGB-DNA on-off processes are fast compared 
with this timescale for both B-form and overstretched 
DNA. In contrast, faster stretch-release cycles reveal 
that the overstretching rate becomes limited by HMGB 
protein dissociation, allowing calculation of the micro- 
scopic dissociation rate. Together with our earlier work, 
these results suggest that HMGB-induced DNA bends do 
not have rigid structures, as HMGB proteins are not 
dissociated by high force. Instead, HMGB binding 
causes site-specific 'softening' of the DNA duplex, and 
the average bending angle characterizes more the degree 
of protein-induced DNA deformation than a particular 
HMGB-bound DNA conformation. Binding events are 
sequence non-specific and transient, with rapid kinetics 
(microscopic residence times ~0.1 10 s) relative to DNA 
looping and other processes where DNA rigidity is 
relevant. HMGB-induced DNA bending is unrelated to 
HMGB-DNA binding cooperativity, which is low. 
Finally, HMGB proteins bind both DNA strands in 
many conformations (including overstretched DNA) and 
inhibit changes in twist on the short timescale of protein 
dissociation. These important features of HMGB architec- 
tural protein function are uniquely revealed through 
the single-molecule optical tweezers approaches described 
here. 
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